# NOT RUN {
season<-c(1,1,2,2,3,3)
team1<-c(1:6)
team2<-c(1,2,5,6,9,11)
team3<-c(10,5,3,7,1,3)
team<-cbind.data.frame(season,team1,team2,team3)
r = c(t(as.matrix(team[,2:4])));
f1 = c("team1", "team2")
f2 = c("season1", "season2","season3");
k1 = length(f1);
k2 = length(f2);
n = 3;
FactD(r,k1,k2,f1,f2,n)
# }
Run the code above in your browser using DataLab